Analysis of major failures in Europe's power grid 
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Power grids are prone to failure. Time series of reliability measures such as total power loss 
or energy not supplied can give significant account of the underlying dynamical behavior of these 
systems, specially when the resulting probability distributions present remarkable features such as 
an algebraic tail, for example. In this paper, seven years (from 2002 to 2008) of Europe's transport 
of electricity network failure events have been analyzed and the best fit for this empirical data 
probability distribution is presented. With the actual span of available data and although there 
exists a moderate support for the power law model, the relatively small amount of events contained 
in the function's tail suggests that other causal factors might be significantly ruling the system's 
dynamics. 

PACS numbers: 84.70.+p, 89.75.Hc, 95.75.Wx 



I. INTRODUCTION 

There has been in recent years an increasing aware- 
ness about infrastructure networks security and reliabil- 
ity [U, H, 3- Modern society's functional capacity relays 
on an optimal operation of infrastructure and informa- 
tion networks such as roads, railways, gas and oil pipes or 
Internet. Particularly vital, and at the same time quite 
prone to failure, are electric power transmission networks. 
These are extremely complex engineered systems, com- 
posed of multiple and interconnected elements, whose re- 
liability depends both on each component's behavior and, 
at the same time, on the many different dynamical inter- 
actions that span over and rule the overall connectivity 
of the system. 

Although it is not always the malfunction of 

a power transmission system shows usually itself as a 
blackout. This is a direct consequence of a cascading 
failure involving several of its composing and linking ele- 
ments. This fact turns the study of the details of failures 
in power transmission networks from a traditional engi- 
neering point of view a hard task, if not an impossible 
one most of the times. In order to reduce the inher- 
ent complexity of this detailed approach, some new ways 
have been proposed in recent years. One of them is that 
of ignoring the details of particular failures and to focus 
on the study of global behaviors and dynamics of time 
series with approximate global models. Concepts such 
as criticality and self-organization have been applied to 
characterize blackout data, suggesting that the frequency 
of large blackouts is governed by non trivial distribu- 
tion functions such as power laws and, consequently, that 
power systems are designed and operated near a critical 
point. (For a comprehensive review on this approach, see 
Ref. and references therein). 

This paper analyses for the first time, and as far as we 
know, the statistics of major electric transmission net- 



work events in the European power grid from this afore- 
mentioned complex systems approach. Following essen- 
tially the statistical analysis presented in Ref . , we esti- 
mate the basic parameters of the power-law model, then 
calculate the goodness-of-fit between the data and the 
power law and finally we compare the power law with 
alternative hypotheses via a likelihood ratio test. The 
paper is organized as follows. In section II European 
major events data is presented and explained. In sec- 
tion III blackout data is analyzed. Finally, section IV 
summarizes our main results. 



II. UCTE MAJOR EVENTS DATA 

European power network reliability data can be found 
in the Union for the Co-ordination of Transmission of 
Electricity (UCTE) web page, publicly available from 
2002 onwards in monthly statistics format [6]. The 
UCTE is the association of Transmission System Oper- 
ators (TSOs) in continental Europe and manages data 
from 24 different European countries. Due to the com- 
plexity of events, sometimes involving more than one 
TSO, types of interruptions in the network and short time 
given to gather this information, UCTE major events 
data is somehow limited in its scope and does not provide 
a fully detailed account of some events. It is, nonethe- 
less, the best documented source that has been found. For 
each major event, it summarizes the following informa- 
tion: 

• Country. 

• Substations involved. 

• Reason (R). Broadly classified into four groups: 
(1) overloads (also calculated brakes), (2) failures 
(false operation, failure in protection device or 
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FIG. 1: Cumulative distribution functions P{x) and their maximum likelihood power law fits for the UCTE reliability measures 
energy not supplied, total loss of power and restoration time. 



other element), (3) external (outside impacts and 
very exceptional weather and natural conditions) 
and (4) other or unknown reasons. 

• Energy Not Supplied (ENS). Measured in 
MWh, as loss of energy from the consumption side. 

• Total Loss of Power (TLP). Measured in MW, 
as loss of production from the generation side. 

• Restoration Time (RT). Measured in minutes. 
Note that since ENS and TLP are measured from 
different sides, RT can not be assumed as the ratio 
of ENS over TLP. It can be considered, therefore, 
an independent reliabihty measure. 

• Equivalent Time of Interruption (ETI). De- 
fined as the duration of an interruption in minutes 
multiplied by the energy not supplied divided by 
the consumption for the last 12 months. Defined 
in this way, the ETI allows a direct comparison be- 
tween TSOs in terms of interruptions that occurred 
during a year. 

From 2002 to 2008, both years inclusive, 908 major 
events have been noticed. Due to the complexity of 
events some entries have zero value in one or more of 
their categories. While these zeroed values have not been 
considered, the rest of numerical values are effective mea- 
sures of major events occurred in the UCTE power grid 
and, consequently, they have all been included in order 
to develop the analysis presented in this paper. 

III. PROBABILITY DISTRIBUTION ANALYSIS 

The study of the statistics and dynamics of series of 
events with approximate global models has been one of 
the most popular topics in the last twenty years, spe- 
cially within the interdisciplinary study of complex sys- 
tems. Probability distribution functions with a heavy 



tailed dependence in terms of event or object sizes seem 
to be ubiquitous statistical features of self-organized nat- 
ural and social complex systems |7j. The appearance 
of algebraic distributions, specially power laws, is often 
thought to be the signature of hierarchy, robustness, crit- 
icality and universal subjacent mechanisms Electric 
power transmission networks have not escaped this cap- 
tivation for power laws quest. Time series of usual mea- 
sures of blackout size like energy unserved, power loss or 
number of customers affected, have been shown to be al- 
gebraically^distributed in North America Sweden p^ . 
Norway ll]. New Zealand [12] and China This ap- 
parent ubiquitous evidence have led to believe and try to 
demonstrate that power systems (a) tend to self-organize 
near a critical point and (b) that there may be some uni- 
versality ruling the inner depths of these systems. 

In spite of this evidence, most of the aforementioned 
literature relay on poorly performed statistical analysis 
and results can not be trusted. In some cases method- 
ologies are not clearly explained (i.e., Ref.[llj) or sim- 
ple visual inspection can clearly dismiss the analysis per- 
formed to rule in the power law hypothesis (i.e., Ref.fToj 
and Ref.[l3]). In other cases, proper usage of statistical 
tools have given new results that limit the scope of the 
original analysis (i.e., Ref.Q is dismissed as insufhciently 
substantiated in Ref. [l^ and reanalyzed in Ref. [1| , find- 
ing moderate support for the power law hypothesis and 
even some for an exponential distribution). 

In this section we analyze the probability distributions 
of three malfunction measures of the European power 
grid: energy not supplied (ENS), total loss of power 
(TLP) and restoration time (RT). The results are sum- 
marized in TableUand shown in Figure[TJ The methodol- 
ogy that has been used is that described in Ref. [5] , where 
a maximum likelihood approach is proposed to estimate 
the heavy tailed function from the data and a signifi- 
cance test is constructed for testing the plausibility of 
the power law. Measures shown in Table [T] are generic 
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a 


ntail 


P 


Support for PL 


ENS 


690 


552 


7004 


180000 


185±72 


1.7±0.1 


104±120 


0.24 


Moderate 


TLP 


576 


400 


1790 


26746 


615±244 


2.1±0.2 


57±96 


0.36 


Moderate 


RT 


897 


510 


3328 


44640 


150±68 


1.69±0.07 


157±115 


0.73 


Ok 



TABLE I; UCTE major events generic statistics and power law fits. For each measure we give the number of occurrences n, 
mean {x), standard deviation a, maximum observed occurrence Xmax, lower bound to the power law behavior Xmin, scaling 
parameter value d, occurrences in the power law tail ntaii and p value p. The last column indicates the support for whether the 
observed data is well approximated by a power-law distribution. Estimated uncertainties for Xmin, & and ritaii are also shown. 



statistics on one side and results of the aforementioned 
statistical analysis on the other. We assume a quantity 
X follows a power law if it is drawn from a probability 
distribution p{x) oc a;~", where a is the scaling param- 
eter of the distribution. Since the probability density of 
a power law distribution diverges as x —^ 0, there must 
exist a lower bound to the power law behavior p^ . We 
denote this lower bound as Xmin and the number of events 
contained in the upper range as ntaii ■ Finally, the p- value 
denotes the significance test result: the power law is ruled 
out if p < 0.1. As we can see, the power law model is 
a plausible one for every data set considered (i.e., the 
p- value for the best fit is sufficiently large) and the scal- 
ing parameter values are similar to those encountered in 
the literature for the ENS and TLP distributions 
Yet the power law model explains only a small amount 
of events: 15% for ENS {ritaii = 104), less than 10% for 
TLP {ntail = 57) and 17% for RT {nta^l = 157) (even 
though it holds the better fit, with p = 0.73). We believe 
that measures such as ntaii and Xmin are fundamental to 
estimate the span of the power law behavior and to de- 
velop further quantitative models, yet these values have 
not been considered in any of the aforementioned refer- 
ences. Only in the reanalysis of Ref.[9i] done in Ref.j^ we 
have found an estimate for ntaii that gives an explanation 
for a 28% of the events. 

We assume that the limited span of available data in 
each set might have a sensible influence in the final power 
law fitting outcome. It is nonetheless evident from these 
results that (a) power law behavior can not be assumed 
for the whole data observed, (b) we can not accept the 
existence of any critical point at this stage of the data 
span and (c) there must be right now considerably more 
dynamics not explained by the power law model. 

In order to check if other distributions may be a bet- 
ter fit, we have performed log likelihood and value tests 
with respect to log-normal, exponential, stretched expo- 
nential (WeibuU) and power law with cut off distribu- 
tions. Results are shown in Table [Til Positive log likeli- 
hood values favor the power law hypothesis and p- values 
higher than 0.1 imply no significance on the results. Log- 
normal, WeibuU and exponential distributions can be 
ruled out as p > 0.1 in the first two and p = in the 
latter with positive LR values. Power law with cut off 
can be ruled out in the ENS and TLP data sets though 



it is a plausible option for the RT data set. 



IV. SUMMARY AND DISCUSSION 

Power outages are considered unexpected phenomena 
in power grids. They appear without warning and, 
though widely investigated, there is not a common ac- 
cepted theory that explains neither their pervasiveness 
nor their inner dynamics. The statistical overabundance 
of big blackouts has been explained using theories of sys- 
tems failure able to reproduce their empirically found 
probability distribution. This distribution is considered a 
power law for most of the literature encountered, with the 
consequences that this algebraic tail involves (i.e., self- 
organization, criticality and universalities). In order to 
add one more reference to this field, in this paper we have 
analyzed seven years of disturbances data for the UCTE 
power grid and for three major event measures: energy 
not supplied, total loss of power and restoration time. Al- 
though evidences for self-organized criticality have been 
suggested for even five years of data [lH| , support for the 
power law hypothesis has been found moderate for two 
(ENS and TLP) of the three measures considered. More- 
over, the amount of events explained by the power law 
hypothesis can be considered negligible. This facts make 
it difficult to accept the existence of an equilibrium point 
near criticality for the UCTE power grid, at least at this 
stage of data analysis, and it also suggests that most of 
the power grid dynamics should be explained by different 
models, other than power law. 

There still exist many complexities not explained in 
this systems. Ongoing research is now focused in analyz- 
ing major events probability distributions in connection 
with (a) the reasons that trigger these major events and 
(b) the structure and topology of the power grids involved 
in these major malfunctions. Results will be published 
elsewhere. 
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log-normal 
LR p 


exponential 
LR p 


Weibull 
LR p 


PL + cut ofr 
LR p 


Data set 


ENS 
TLP 
RT 


-0.405 0.68 
0.319 0.75 
-0.382 0.70 


2.64 0.00 
3.42 0.00 
7.57 0.00 


-0.393 0.69 
0.467 0.64 
-0.329 0.74 


-0.419 0.36 
-0.08 0.68 
-1.7 0.06 



TABLE IL Test of power-law behavior. Positive values of the log likelihood ratios LR favors the power law model. Values of 
p > 0.1 imply though that result can not be trusted. The exponential distribution is definitely ruled out as possible model and 
only for the restoration time RT, the power law with cut off could be considered a valid model. 
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